LANDAU-95-GR-1 



^1" 



MODERN COSMOLOGICAL MODELS WITH DARK MATTER 
AND THEIR CONFRONTATION WITH OBSERVATIONAL DATA 

A. A. Starobinsky 
Landau Institute for theoretical Physics, 
Russian Academy of Sciences, Moscow 117334, Russia 

. Short version of the plenary talk presented at the First International Conference on Cos- 
moparticle Physics "Cosmion-94" (Moscow, 5-14 December 1994)- Published in: Cos- 
moparticle Physics. 1, eds. M.Yu. Khlopov, M.E. Prokhorov, A. A. Starobinsky, J. Tran 
Thanh Van, Edition Frontiers, 1996, pp. 141-146. 

SI 

ABSTRACT. New systematic classification of cosmological models of the present Uni- 
verse is introduced. After making the comparison of these models with all existing obser- 
vational data three viable models remain: the cold dark matter model with the cosmo- 
logical constant (which becomes the most reasonable one if the Hubble constant Hq > 60 
km/s/Mpc); the mixed cold-hot dark matter model (models with two and especially three 
types of neutrinos with equal masses are in a slightly better agreement with observational 
data than the model with one massive neutrino); the pure cold dark matter model with 
a step-like initial spectrum of perturbations. The two latter models require H < 60 
km/s/Mpc. 
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^p" 1 . Looking at the enormous amount of papers in which different cosmological models 

of the present Universe are studied and compared with observational data, the first im- 
pression might be that there exists a great degree of freedom and arbitrariness in such 
models. However, actually this is not so. First, using a natural fundamental classification 
described below, it appears that all currently popular models follow some simple logic. 
It becomes clear also what are the natural directions of generalization of these models. 
Second, if all existing observational data are used (and not some arbitrarily chosen tests 
only), then a number of viable cosmological models reduces to a few - we shall argue below 
that only three models remain viable at present. Would we know the Hubble constant Hq 
with a better accuracy, we could further reduce this number and might even finish with 
a unique preferred one. 

To introduce this classification, let us remind the logical structure of cosmology. The 
aim of theoretical cosmology, as of any other branch of science, is to produce definite 
predictions basing upon as minimum number of assumptions as possible (" Occam's razor" ) 
which can be tested by observations. To make these predictions, first, we need equations 
(or Lagrangians) . It is generally accepted now that the Universe is described by a space- 
time metric satisfying the Einstein equations with some matter source. Moreover, it was 
shown that a broad class of theories more general than the Einstein gravity (e.g., scalar- 
tensor theories of gravity, the R + f{R) theory where R is the trace of the Ricci tensor, 
etc.) can be represented in the Einsteinian form after some conformal transformation. 
So, the left-hand side of the equations may be taken as Rik — \Rgik- 

However, the equations are not defined if their right-hand side is not completely speci- 
fied. So, we should assume some matter content of the Universe (in particular, at present 
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time). It is remarkable that all viable cosmological models require the main part of matter 
in the present Universe to be the dark non-baryonic matter. It is dark in the direct sence 
of the word, i.e. not participating in electromagnetic interactions. Moreover, all candi- 
dates for this dark matter are not yet discovered in laboratory experiments (or not yet 
discovered are the specific properties of particles making them important for cosmology, 
if speaking, e.g., about neutrino masses). Thus, modern cosmology strongly suggests and 
awaits great laboratory discoveries in the field of particle physics! 
The evidence for dark matter follows from the following numbers. 

1. The energy density of luminous matter in the Universe (stars, hot gas in galaxies 
and between them) in terms of the critical one e c = 3Hq/8tiG (we assume c = H = 1): 
Vt lum < 0.01. 

2. The primordial (Big Bang) nucleosynthesis (BBN) prediction for the energy density 
of baryons in the Universe based on observed abundances of light elements ( 4 He, D, 3 He, 
7 Li) is tt b = (0.01 - 0.02)h- 2 where h = H o /100. 

3. The virial energy density of matter gravitationally clustered at scales R = (5 — 25) /i -1 
Mpc is tt vir . = 0.2 - 0.5. 

4. Finally, the theoretical expectation in most inflationary models is fl t ot = 1- 
Comparing 2 with 1, we see that some amount of baryon dark matter should exist (e.g., 
in the form of faint stars). However, from the comparison of 3 to 1 we see that much more 
matter should be in a non-baryonic dark form. Note that this conclusion is not based on 
any theoretical assumptions about a cosmological model. If we compare 4 with 3 (using 
now some hypothesis about the early Universe), then the further conclusion follows that 
there may exist different types of dark non-baryonic matter in the Universe. 

The main 3 types of dark matter (further, we shall speak about non-baryonic dark 
matter only) are: 

1. Hot dark matter (HDM) - neutrinos with restmasses of the order of few eV. 

2. Cold dark matter (CDM) - sypersymmetric particles with masses ~ 100 Gev or axions 
with m a ~ 10~ 5 eV. 

3. Ultra-cold (vacuum-like) matter - the cosmological constant. 

However, we understand now that there are no impenetrable barriers between these 3 
types. In more exotic models, intermediate types of dark matter may appear between 
these ones. In particular, warm dark matter (e.g., thermal particles with masses of about 
1 keV or non-equilibrium neutrinos) fills the gap between 1 and 2, and a scalar field with 
the exponential potential can be used to produce an arbitrary equation of state of matter 
between 2 and 3 (actually, even in the range from p = e to p = — e where p is the pressure 
of matter). 

So, if the dark matter content is specified, the equations are fully determined. However, 
to solve them we need initial conditions (mostly, for perturbations). The observational 
evidence from AT/T angular anisotropy of the cosmic microwave background (CMB) 
proves that the Universe was FRW-like (isotropic and homogeneous) at the recombination 
time (at the redshift z ~ 1100). The BBN theory shows that the Universe was FRW-like 
in the much earlier period (t = (1 — 100) s after the cosmological singularity, z = 10 8 — 10 9 ). 
The absence of primordial black holes with M = 10 15 — 10 17 g which would evaporate 
now through the Hawking radiation suggests the isotropy of the Universe at even earlier 
times, and so on. However, even if we assume that the Universe was FRW-like just from 
the very beginning, i.e. from times close to the Planck time tp = \/G, we still remain 
with three arbitrary functions of space coordinates (or 3 arbitrary functions of the wave 
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vector k in the Fourier representation) specifying initial amplitudes of the quasi-isotropic 
modes (called so because they do not destroy local isotropy and homogeneity at arbitrarily 
early times). One of these functions refers to scalar (adiabatic) perturbations and two 
others - to the quasi-isotropic mode of tensor perturbations (primordial gravitational 
waves). In more complicated models, additional functions giving initial values for non- 
decreasing isocurvature modes can appear. Therefore, generally we have a large functional 
arbitrariness in initial conditions. 

The second main advantage of the inflationary scenario of the early Universe, after 
the elegance and beauty of its main assumption that our Universe was in the maximally 
symmetric (de Sitter, or inflationary) state during some period in the past, is that it 
predicts these initial conditions in terms of a small number of fundamental parameters of 
an effective Lagrangian describing an inflationary stage realized in any concrete version of 
this scenario. Thus, the arbitrariness reduces to a few (minimum one) parameters. From 
all this discussion, a natural and fundamental classification of cosmological models of the 
present Universe follows (including non-inflationary models, too, but I shall discuss only 
inflationary ones further) [1] (see also [2]): 

let us classify them by their level of complexity number which is equal to a number of 

a) new (not known before, e.g., H and f2 7 are not counted), 

b) fundamental (appearing in basic equations, not in initial conditions), 

c) significant (more than ~ 10% effect, e.g. if \n s — 1| < 0.1, it is counted as n s = 1), 

d) dimensionless (this can be always achieved, e.g., working in the Planck units) 
constants introduced in any particular models in order to explain all observational data. 
These constants may refer either to the present dark matter content or to the initial 
spectrum of perturbations. Note that the complexity level of a given model may decrease 
as a result of future discoveries in laboratories (e.g., measurements of neutrino restmasses) 
or due to a progress in an underlying unified physical theory (e.g., a derivation of inflaton 
parameters from the superstring theory). 

Distinguishing features of this classification which radically discriminate it from all 
previous attempts in this direction are, first, the fact that the number of assumptions 
using no numbers (numerical constants) is not counted at all. Thus, this classification is 
favourable even for very "strange", "crazy" models if they are internally consistent and do 
not introduce a large number of free parameters. Second, it counts assumptions referring 
to the present dark matter content and to the initial spectrum of perturbations on equal 
footing. Thus, if we want to complicate a model under the pressure of observational data 
and shift to the next level of complexity, we may add one parameter eigher to the effective 
Lagrangian describing an inflationary stage (and then we get one more parameter in the 
initial spectrum of perturbations), or to the description of the dark matter content (e.g., 
by introducing one new type of dark matter). Then the logical order of the development 
is to begin with the lowest complexity level, compare it with the data, and if there is no 
agreement, move to successive higher levels until the agreement will be reached. 

In the Table 1 below a very brief sketch of existing cosmological tests is presented 
with an approximate range of scales to which they are sensitive (or of scales for which 
data exist at present). Here $ is the gravitational potential and h a p are the tensor per- 
turbations (gravitational waves). The Table 2 gives a list of models with their complexity 
level numbers, initial (fundamental) parametres, observational parameters which can be 
unambiguously expressed through the fundamental ones, and the result of comparison 
with observational data. A is the rms amplitude of the spectrum of adiabatic perturba- 
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tions, Ai SO is the same for isocurvature perturbations, <jg is the total rms top-hat matter 
perturbation at the scale R = 8/i _1 Mpc. For the model 6, 0o is the value of the inflaton 
field at the moment of bubble formation, 0/ is its value at the end of inflation. For the 
model 10, <p s is the value of the inflaton field at the beginning of the last stage of inflation. 
For the model 11, v\ denotes an unstable neutrino while v denotes the stable one. The 
latter model (which is actually a class of different models) is still largely unexplored, so 
no conclusion about its viability is drawn. 

All details about the confrontation of the models with observational data can be found 
in [3,2] for the models 3 and 9, in [4,5] for the model 4, in [6] for the model 5 and in [7,8] 
for the model 10. The list of references is, of course, very incomplete, but references to 
other papers can be found in the given ones. 

So, the first remarkable conclusion is that it is possible to explain all existing cosmo- 
logical observational data using 2 or 3 parameters only (4 in the case of three neutrinos 
with comparable masses). There no necessity to go to higher complexity levels at present. 

The second conclusion shows the dependence of results on the value of H Q (still not 
determined with the desired degree of accuracy). If H > 60 km/s/Mpc, then the best 
(and probably the only possible) model is the model 5: CDM+A with n s ~ 1. So, a 
reliable observational proof that Hq > 60 should be considered as a very strong argument 
for the positive cosmological constant. Note, however, another prediction that Hq < 80 
because in the opposite case no reasonable model exists (at least, at the complexity levels 
considered). On the other hand, for H Q < 60 we can avoid introducing the cosmological 
constant, and then the choice is between the mixed CDM+HDM models 3 or 4 and the 
pure CDM model 10. They can be most easily discriminated by the abundance of compact 
objects at large redshifts. 

The third conclusion is that if the mixed CDM+HDM model is the right one, then it is 
very interesting that the model 4 having several neutrino species with comparable masses 
produces a fit to all data slightly better than the model 3 with one massive meutrino. So, 
in this case cosmology provides some support (though not the final evidence, of course) 
for strong neutrino mixing. 

This research was partially supported by the Russian research project "Cosmomicro- 
physics" through Cosmion and by the INTAS grant 93-3364. 
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Table 1. Tests 



Name 


Type of perturbation 


Scales (h 1 Mpc) 


Large-angle AT/T (6 > 2 U ) 


&,h afi 


200-6000 


Intermediate angle AT/T (10' < 9 < 2°) 


$, V$, A$ 


20-200 


Peculiar (bulk) velocities 


V$ 


10-75 


Galaxy-galaxy correlations 


A$ 


0.1-200 


Large-scale structure 


a 2 * 


5-100 


Cluster-cluster correlations 


A$ 


5-150 


Cluster abundance 


A$ 


5-10 


Galaxy and quasar abundance at large z 


A$ 


0.5-1 


Ly — a clouds 


A$ 


0.1-1 
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Table 2. Classification of cosmological models 



No. 


Models 


Level No. 


Fundamental 
parameters 


Observational 
parameters 


Agreement with 
observational 
data 


1 


Standard CDM 


1 


(H 2 /<p) inf 


A 


No 




(n s = 1) 








2 


CDM, n s = -3 
isocurvature 


1 




/t . 


No 


3 


CDM+HDM, 

lZA 71, = 1, 

7 £> 7 

standard 
concentration 


2 


(# 2 /0W, 
m„lM n 

Is 1 P 


A, Vt v 


Marginally 
good if 
H ~ 50, 
Q„ ~ 0.2 


4 


CDM+HDM 
2z/ or 3z/, 
n s = 1, 
standard 

concentration 


3,4 


(H 2 /6) inf , 

\ It) iilj i 

Tfly-^ 771 

(= m u 3 ) 


A.* 


Good if 
iJ < 60, 
Em, ~5(ifo/50) 3 eV 


5 


CDM+A, 

Qtot — 1, 
n s = 1 


2 


(H 2 /<P) mf , 
AM 2 


A,n A 


Good, 
H = 50 - 80, 
fi m = 0.5 - 0.2 


6 


CDM curved, 
n s = 1 


2 


(H 2 /<j>) inf , 

00 l<t>} 




Worse than 
CDM+A, 
n m > 0.3, 
H < 70 


7 


CDM 


2 


V(d>) = Vne~ a ^, 


/■ 1 it lit s\v 1 - 

- 1 x \ ,K TlOV J 1 


No 




tilted 




cy Vn 

, r (J 


TI,o ^ 1 


(if cr« < 7) 




adiabatic 










8 


CDM 


2 






No 




tilted 






n s 7^ -3 






isocurvature 










9 


CDM+HDM 


3 


a, Vb, 


A(hhor), 


No, if 




tilted, lv, 




m u /M p 


n s ^ 1, 


\n s - lj > 0.1 




standard 






n„ 






concentration 










10 


CDM with 


3 


(H 2 /<j>) + , 


A + ,A_, 


Good if 




a step-like 






k s 


# < 60 




spectrum 




<f>s/</>f 

(H 2 /<P) mf , 






11 


CDM+HDM, 
decaying u, 
n s = 1 


3 


A, 
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